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ABSTRACT 

‘iwo  He-Ne  lasers  operating  at  6328  %.  have  been  utilized  to  make 
simultaneous  measurements  of  the  effects  of  scintillation  over  homogeneous 
optical  paths  of  650  and  1300  m  to  study  the  transfer  of  coherent  radi¬ 
ation  through  a  turbulent  medium.  At  the  path  terminus,  multiple  sampling 
of  each  laser  beam  was  effected  using  a  photo-optical  technique  which 
records  a  61-cm  cross  section  of  an  optical  beam.  Concurrent  with  the 
optical  data,  wind-speed  and  direction  recordings  were  made  at  multiple 
points  along  the  optical  path  in  order  to  estimate  the  homogeneity  of 
meteorological  conditions.  Near  the  path  terminus,  measurements  of  wind 
shear  and  temperature  lapse  were  taken.  In  addition,  high-speed  ther¬ 
mometry  techniques  were  utilized  to  compute  one-dimensional  temperature 
spectra  as  well  as  the  thermal  structure  coefficient  Cj.  Data  were 
gathered  during  temperature-lapse,  neutral,  and  inversion  condition'*. 
Log-irradiance  scans  derived  from  the  optical  data  were  used  to  compute 
log-irradiance  power  spectra,  variance,  and  other  statistical  quantities. 
Using  these  optical  and  meteorological  data,  optical  filter  functions 
were  then  calculated  for  spatial  frequencies  above  87  c/m  and  are  used 
to  compare  with  current  theories.  The  saturation  of  the  log-irradiance 
data  is  again  observed,  and  the  isotropy  of  the  irradiance  fluctuations 
is  examined. 
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INTRODUCTION 


Light  propagation  through  the  lwer  atmosphere  has  r«c«iv«d  eon- 

siderable  attention  in  the  past  decade,  both  theoretically  and  sxperi- 

1* 

mentally.  A  recent  review  article  by  Strohbehn  listt  5?  references 
on  the  subject.  Therefore,  a  detailed  historical  sisatary  of  the  sub- 
ject  is  omitted,  but  selections  from  current  theoretical  developments 
directly  relevant  to  this  paper-  are  included. 

It  is  well  known  that  a  light  wave  propagating  through  •  region 
of  turbulent  atr-oopiiave  suffers  phase  and  amplitude  distortions  as 
it  progresses.  If,  instead  of  an  infinite  plane  wave,  a  narrow  been 
is  propagated  through  the  medium,  the  center  of  energy  of  the  bea* 
can  be  observed  to  wander  and  the  beam  diameter  to  vary  iu  tine..  The 
beam  is  being  steered  by  passing  through  a  set  of  spatially  and  tem¬ 
porally  random  refractive- index  gradients.  If  the  diameter  of  the 
beam  is  increased,  marked  irradiance  variations  (scintillations) 
become  evident  at  sufficient  ranges  from  the  transmitter,  while  the 
beam  steering  effects  are  gradually  reduced.  For  certain  path- 
turbulence  conditions,  spectacular  optical  effects  can  occur,  yielding 
large,  irregular  beam  patterns.  Finally,  as  the  transmitted  beam 
assumes  a  spherical  wave  front,  the  scintillation  becomes  independent 
of  the  transmitter  and  a  function  only  of  the  spectral  character  of 
the  atmospheric  turbulence. 

References  are  listed  on  page  42. 
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The  three  classes  of  optical  effects — beam  wander,  diameter 
variation,  and  scintillation — are  related  directly  to  the  low, 
medium,  and  hign  wave  number  regions,  respectively,  of  the  turbu¬ 
lence  spectrum.  The  relative  magnitudes  of  these  three  classes  of 
optical  effects  are  closely  related  to  the  transmitter-range  parameters 
as  well  as  to  the  variant  character  of  the  turbulence. 

THEORETICAL  BACKGROUND 

Current  light  fluctuation  theories  may  be  classed  according  to 

whether  single  or  multiple  scattering  is  considered  and  whether  the 

initial  unperturbed  wave  is  eithei  an  infinite  plane  (or  spherical) 

wave  or  that  of  a  finite  beam.  Of  the  single-scattering  theories, 

2 

the  one  formulated  by  Tatarski  is  perhaps  the  most  widely  quoted. 

Tatarski  considers  a  plane  or  spherical  wave  impinging  on  a  slab  of 

turbulent  atmosphere  whose  statistics  are  known.  For  many  systems 

applications,  an  actual  source  might  be  better  described  as  a 

diffraction- limited  beam  with  a  gaussian  transverse  irradiance  pro- 

3  4 

file.  Theories  put  forth  by  Ishimaru,  Carlson,  and  Tatarski  and 
Kon^  describe  such  beam  sources. 

As  might,  be  expected,  the  introduction  of  an  additional  finite 
transmitter  aperture  changes  the  predicted  statistics  of  the  received 
light.  Scintillation  is  produced  by  interference  between  scattered 
and  unscattered  portions  of  the  beam.  If,  for  a  given  range,  the 
beam  is  sufficiently  narrow  so  that  the  phase  shifts  between  beam 


elements  are  limited  to  much  less  than  n  radians,  then  irradiance 
variations  may  be  largely  eliminated.  The  predictions  of  phase  and 
scintillation  spectra,  as  well  as  phase  and  log-irradiance  variances 
calculated  by  plane-wave  single-scattering  theory  (PST)  and  beam-wave 
single-scattering  theory  (BST),  are  subject  to  experimental  test  as 
will  be  described  below. 

If  the  standard  error  of  the  measured  log-irradiance  of  scintil¬ 
lation  is  plotted  against  the  theoretical  prediction  obtained  from  PST 
for  given  turbulence-range  parameters,  it  has  been  shown  by  Gracheva 
and  Gurvich^  in  the  Soviet  Union,  and  verified  by  Deitz  and  Wright ^ 
and  others  in  this  country,  that  the  statistics  "saturate."  That  is, 
the  measured  statistics  become  nearly  independent  of  the  turbulence- 
range  parameters  beyond  a  certain  critical  value  of  the  latter.  Expla¬ 
nation  of  this  phenomenon  has  stirred  interest  in  multiple-scattering 
calculations,  the  chief  contributors  being  Tatarski  and  deWolf. ^ 

The  plane-wave  multiple-scattering  theories  tPMSl)  cited  above  do  give 
rise  to  saturation  curves;  the  one  most  cited  is  due  to  deWolf.  How¬ 
ever,  the  complex  nature  of  these  formulations  has  limited  to  date 
their  application  to  predictions  of  log-irradiance  variance;  the  trans¬ 
verse  phase  power  spectra,  for  example,  are  not  given  by  current  PMST. 

It  is  important  to  note  that  if  the  variance  predicted  by  BST  is 

plotted  against  the  results  predicted  by  PST,  a  linear  relationship  is 

not  revealed  either.  The  general  features  of  the  saturation  measure- 

4 

ments  seem  to  be  described  by  the  formulation  of  Carlson,  but  with 
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rather  significant  departures  predicted  at  long  ranges.  Current 
experimental  work  is  not  yet  adequate  to  make  distinctions  at  long 
ranges . 

We  shall  confine  our  comparison  of  the  several  theories  to 

(3) 

examination  of  scintillation.  Let  <i>^  (x)  be  the  isotropic  three- 

dimensional  power  spectrum  of  the  atmospheric  refractive-index  fluc¬ 
tuations  at  the  spatial  wave  number  k-  Let  k  be  the  optical  wave 
number,  L  the  path  length,  0(x)  the  optical  filter  function  of  the 
atmosphere,  x  the  (random)  log-amplitude  of  the  light  wave,  a  -  A/irr0 

the  inverse  Fresnel  distance  of  an  aperture  of  radius  r  ,  and  F  (x,0) 

X 

the  two-dimensional  power  spectrum  of  the  beam  cross  section 
referred  to  object  space.  These  quantities  are  related  by 


F  (k,0)  =  7rk2L  9(<)  4>^(x) 

A 


(1) 


for  all  first-order  scattering  theories.  Now  it  is  also  true  that 
2 

the  variance  <x  >(<x>=  0)  is  given  as  the  integral  of  F  where 

X 

CO 


<X2>  =  2  it 


F  (x,0)  k  dx. 


(2) 


The  variance  may  also  be  obtained  by  more  direct  averaging  as  in  the 
case  of  PMST. 

2  (3) 

The  quantities  <x  >F  (k,6)  and  4>  (x)  are  all  accessible  through 

X  n 

measurement.  Consequently,  it  is  possible  to  obtain  the  optical  filter 

function  by  means  of  Eq.  (1)  and  to  compare  this  result  along  with  the 
2 

measured  <x  >with  those  predicted  by  theory.  The  formulations  are 

12 

given  below.  The  PST  gives 
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0(0  =  i  sin[  k  ]  ’ 

(3) 

2  ,  -2  .  7/6  , 11/6 

<  >  "  0.307  Cn  k  L  , 

(4) 

3  4 

whereas  BST  ’  predicts 


The  variance  in  the  BST  is  given  by 
<x2>  =  7,2  0.033  C2  r (-5/6)  k7/e  Lll/6[i-r^L?]5/6[|-  *(«L>]  ,  (7) 

where 

g(c,L)  =  Ke  [jf  ,F,(=|  i  1  :  ^  i  id)]  •  (8) 

Finally,  FMST  predicts  a  variance  of 

<X2>=  \  In  |^2  -  exp(-4  <y.^>)j.  v9) 

Plots  of  the  optical  filter  functions  described  by  Eqs.  (3)  and 
(5)  are  shown  in  comparison  with  experimentally  derived  curves  at  650- 
and  1300-m  ranges  in  Figs.  1  and  2.  The  log-irradiance  statistics 
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RANGE 
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f (CYCLES  PER  METER) 

Optical  Filter  Function,  o(f),  as  a  Function  of  Spatial  Frequency 
predictions  of  Tatarski  (TFF)  and  Carlson  and  Ishimaru  (CFF)  are 

rn  in  solid  lines,  and  experimental  results  for  several  times  are 


RANGE 
1300  M 


v 


>> 

u  QJ  gj 

C  U  U 
Hi  «j  f a 

QJ  U-  CJ 
Ui.  e 
LU  O  *r= 


3  h— 

•r-  nj 
+->  <0  t- 

e  oj 

C-'r-  > 

oo  .c  <u 

,  m  in 

V-  h-« 

o  u 
*D  o 
C  C  4- 
O  ro 
•i-  tn 
P  C+J 
O  Oh-  . 
c  in  3  E 
3  r—  in 
o_  i-  a;  o 
(0  uo 

(OU  CO 
f—  v— 
in  *T3  fQ 
fO  C  M  4- 
IQ  C  O 
*  ,0) 

* — *<-*>  t;  QJ 
H-  U.  •*- 
— 'LU  t-  C 
o  h-  CJ  IT3 

^  C.  U 
*  X 
C  *r  QJ 
O  -X 

•r  JOt  t 

•m  *.  r.  c 

o  nj  <o  H- 

C  4-> 

3  fl  *  in 

u.  p-  m  <u 

<1)  C 
fc-  4-  C  -r- 
(D  O-rr- 
4J  r- 
*—  CO  C 

c-o  <d 

U.  O  »r- 

•r-  r-  O 
•—  4->  O 
nj  U  in  n 
<J  *1— 

•r-  TO  C  C 

4-)  QJ  *r-  «r- 

Q.  t- 

oc.cc 

QJ  QJ  O  O 

r 

K—  in  in 


csi 

QJ 

I- 

3 

cr> 


are  presented  in  Fig.  3.  In  the  latter,  twice  the  square  root  of 
2 

<  XT  >  is  plotted  along  the  abscissa,  and  the  other  theoretical 
results  [Eqs.  (4),  (7),  and  (9)]  are  shown  as  appropriately  identified 
curves.  The  spread  in  the  data  and  the  uncertainty  in  the  measure¬ 
ments  are  shown  by  the  error  curves. 

THE  EXPERIMENT 

Full  documentation  of  an  optical  propagation  experiment  suggests 
detailed  measurement  of  high-frequency  temperature  and  velocity  sta¬ 
tistics  at  multiple  points  along  the  optical  path.  Such  an  approach 
is,  of  coursa,  limited  by  practical  considerations.  The  ideal  optical 
experiment  is  conducted  over  a  uniform  terrain  above  which  statis¬ 
tically  homogeneous  turbulence  is  maintained.  For  such  conditions, 
detailed  measurement  of  the  micro temperature  structure  and  fluctuation 
power  spectra  need  be  made  at  only  one  point.  A  useful  compromise 
under  practical  conditions  is  to  make  the  detailed  measurements  at 
one  or  more  points  near  the  beam  and  to  record  the  gross  meteoro¬ 
logical  conditions  contained  in  wind-speed  and  direction  measure¬ 
ments  at  a  number  of  points  and  altitudes  along  the  path.  The  gross 
measurements  thus  serve  to  test  for  homogeneity  of  conditions. 

The  BRL  Optical  Propagation  Range,  at  which  the  experiments  were 

conducted,  is  a  leveled  grass  field.  Local  obstacles  were  removed  to 

13 

maintain  homogeneous  turbulence  for  specific  wind  directions.  A 
650-m  section  of  the  optical  range  was  instrumented  at  152-m  intervals 
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PLANE  WAVE 


Figure  3.  Standard  Errors  of  Log-Irradiance  as  a  Function  of  the  Plane  Wave  Tatarski  Prediction 


with  wind-speed  and  direction  sensors  and  low-frequency  thermometers. 

At  each  of  five  locations,  a  set  of  sensors  was  mounted  on  poles 
adjacent  to  the  optical  path  at  heights  of  2.5,  5.2,  and  9.4  m. 

The  outputs  of  these  sensors  were  selectively  scanned  and  recorded 
by  an  analog-to-digital  conversion  system.  An  additional  set  of 
five  wind-direction  indicators  having  capability  of  real-time  readout 
was  installed  adjacent  to  the  optical  path  so  that  optimum  times 
could  be  chosen  for  data  acquisition. 

As  evident  by  Eqs.  (4)  and  (7),  the  refractive-index  structure 
constant  Cn  is  the  measure  of  turbulence  power  entering  the  expression 
of  scintillation  variance.  Provided  an  "inertial  subrange"  of  turbu¬ 
lence  exists  described  by  a  Kolmogoroff  velocity  spectrum  and  that 

temperature  inhomogeneities  can  be  treated  as  conservative  passive 
14 

additives,  then  the  mean-square  fluctuation  of  temperature  is  given 
by 

2/3 

<  [T(rj)  -  TCl,)]  >  s  DT(  |  -  rj)  -  Gp  |rj  -  |  .  (10) 

Electrical  measurement  of  the  left  side  of  Eq.  (10)  by  means  of  cold- 
wire  resistance  thermometers  at  two  fixed  points,  ^  and  and 

suitable  difference-squaring  and  averaging  leads  to  the  universal 
thermal  structure  constant  C^,.  For  slow  density  variations  charac¬ 
teristic  of  turbulence,  the  atmosphere  is  approximately  incompressible. 
Therefore,  the  refractive-index  structure  constant  is  given  by  a 
simplified  Barrel-Sears  formula^  where 
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(11) 


where  p  is  the  atmospheric  pressure  in  millibars  and  T  is  the  abso¬ 
lute  temperature. 

The  temperature  structure  constant  may  also  be  inferred  by 
measurement  of  the  vertical  wind  velocity  and  temperature  profiles 


by  the  computation  of  the  approximate  Richardson  number  where 


g(Tl  -  T2)  (Zlz2)' 
Ta„b(?l  -  V 2 


log (s1/z2) ;  g  =  980  cm/sec2' 


An  average  Ri  for  pairs  of  points  measured  is  readily  computed. 
Theory  and  measurements  by  Tsvang^"8  indicate  that 


CT  -  f (Ri) •  (13) 

Values  of  CT  were  derived  through  use  of  Tsvang’s  published  curve. 

Conversions  were  then  made  of  C„  to  C  . 

i  n 

Measurements  of  the  structure  constant  by  both  the  differential- 

probe  and  the  Richardson-number  methods  are  given  in  Table  I.  Probe 

separations  of  10  and  36  cm  were  used  in  the  differential  technique. 

The  vertical  temperature  and  average  wind  velocity  gradients  were 

sampled  at  five  logarithmically  spaced  altitudes  to  a  maximum  height 

of  8  m.  A  second,  and  independent,  measure  of  C^,  as  well  as  the 

17 

spectrum  of  microtemperature  fluctuations,  was  made  utilizing  a 
modified  single  cold-wire  probe.  The  index  of  refraction  spectra, 
derived  from  the  microtemperature  spectra  by  applying  the  Barrel- 
Sears  conversion  factor,  are  shown  in  Figs.  4 A  and  4B. 


19 


METEOROLOGICAL  DATA 
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Figure  4A.  Index  of  refraction  spectra  and  the  corresponding  Ri  estimate 
are  given  above  for  each  curve.  Spectra  were  measured  by 
Portman  and  co-workers.  Associated  times  ara  1130h  (1  and  2) 
and  1430h  (3). 


Micrometeorological  and  optical  data  taken  on  14  June  1968  can 
be  characterized  by  the  vertical  te^oerature  gradient.  The  first  two 
data  runs  (1130  and  1430  h)  took  place  during  temperature  lapse  con¬ 
dition?,  the  third  (1900  h)  during  approximately  neutral,  and  the 
final  two  (2130  and  2200  h)  occurred  during  inversion  conditions. 

The  day  was  clear  with  less  than  five  percent  cloud  cover,  thus  mini¬ 
mizing  variant  solar  absorption  at  the  ground.  The  statistics  for 
the  meteorological  tables  were  computed  from  ten-minute  data  runs. 

The  indicated  times  correspond  to  the  middle  of  the  run  and  the  time 
at  which  the  optical  data  were  taken. 

For  the  optical  measurements,  two  He-Ne  lasers  operating  at 
o 

6328  A  were  utilized.  Each  laser  output  was  collimated  through  a 

telescope  employing  a  spatial  filter  to  remove  variations  in  the  wave 

front.  Eoth  beams  were  then  adjusted  to  a  divergence  of  1.9  mrad, 

since  this  parameter  has  been  found  to  be  a  factor  in  the  resulting 
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irradiance  fluctuations.  ’  One  of  the  lasers  was  located  at  the 
south  end  of  the  propagation  range  and  directed  north  over  the  650-m 
path  to  an  optical  receiver  system.  The  second  laser  was  positioned 
adjacent  to  the  receiver  and  directed  toward  the  south  end  of  the 
range.  An  optical  flat,  located  next  to  the  first  laser,  redirected 
the  beam  northward  to  be  recorded. 

The  optical  receiver  system,  essentially  a  telescope  with  a  61-cm 
primary  element,  is  arranged  with  the  image  plane  conjugate  to  the 
entrance  aperture.  Placing  photographic  film  at  the  image  and  using 
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an  appropriately  short  exposure  time,  a  time-frozen  recording  can  be 
made  of  the  turbulence  fluctuations  in  the  laser  beam.  The  calibra¬ 
tions  and  method  of  data  extraction  by  this  technique  have  been  earlier 
7  18  4 

documented.  *  This  method  has  certain  advantages  over  temporal  sensing 
techniques.  Microdensitometer  scans  of  these  photographs  allow  measure¬ 
ments  to  be  made  directly  in  the  spatial  domain,  thus  avoiding  the 
inherent  approximations  of  converting  temporal  fluctuations  by  use  of 
the  wind  component  normal  to  the  optical  path.  Further,  no  assumption 

19 

need  be  made  about  the  degree  to  which  the  fluctuations  are  "frozen  in." 

In  addition,  scans  may  be  made  through  the  image  in  any  direction, 
including  the  vertical,  making  possible  examination  of  the  isotropy  of 
the  irradiance  fluctuations. 

Five  sets  of  optical  data  were  taken  corresponding  to  the  five 
20 

times  indicated  in  the  meteorological  table.  At  each  period,  three 
successive  photographs  were  made  of  the  laser  beam  cross  section  at 
1300  m  using  an  exposure  time  of  2  msec.  Immediately  following,  three 
photographs,  taken  at  the  same  exposure  time,  were  made  at  the  650-m 
range.  Thus,  the  six  photographs  taken  during  each  run  were  acquired 
over  an  interval  of  about  30  sec.  A  horizontal,  vertical,  and  two 
diagonal  scans  were  made  of  each  photograph  using  a  microdensitometer 
with  an  effective  aperture  (referred  to  object  space)  of  0.4  x  3  mm. 
Corrections  for  this  aperture  are  contained  in  the  composite  measured 
transfer  function  made  for  this  instrument. 
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DATA  REDUCTION 


Digitized  data,  obtained  from  the  densitometer  scans  of  the 
photographs,  were  processed  according  to  the  flow  diagram  given  in 
Fig.  5.  A  step  tablet  with  approximately  thirteen  steps,  photographed 
with  the  optical  data,  established  the  characteristic  curve  of  the 
film.  Previous  measurements  of  the  optical  system  have  indicated  that 
through  the  spatial  frequency  range  of  interest  the  transfer  function 
of  the  optics  is  effectively  unity.  Hence,  the  major  transfer  function 

for  which  to  account  was  that  of  the  densitometer.  The  measured  fune- 

2 

tion  was  fitted  to  a  sine  function  where 

T3(k)  =  sine2 (0.00154k/2)  (14) 


where  <  is  the  spatial  wave  number  in  inverse  meters.  The  standard 
errors  of  the  irradiance  and  the  log-irradiance  data  were  computed  in 
Programs  I  and  II  as  shown  in  Table  II.  Since  for  normally  distributed 
data 

=  ln(l  +  Oj);  1  =  1,  (15) 

2 

where  is  the  variance  of  the  irradiance  data,  it  follows  that  a 
measurement  of  directly  along  with  0^  is  a  measure  of  the  log- 
normality  of  the  statistics.  The  column  in  Table  II  marked  "corrected 
scale  ratio"  gives  the  number  of  irradiance  correlation  intervals 
across  the  scanned  portion  of  the  mirror.  The  defining  expression  is 


Corrected  Scale  Ratio  = 


500  o2 


2 


B:(l)  dl 

0 


(16) 
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Figure  5.  Data  Reduction  Flow  Diagram 

The  individual  subsystem  transfer  functions  are  represented  by 
T  through  T4  and  the  overall  system  transfer  function  by  T. 
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TABLE  II 

STATISTICAL  AVERAGES 


Time 

Range 

°2X  | 

* 

°2X 

Corrected 
Scale  Ratio 

650 

1.036 

1.018 

25.29 

1130 

1300 

0.676 

0.811 

30.35 

650 

1.042 

1.050 

24.82 

1430 

1300 

0.682 

0.823 

31.81 

650 

0.498 

0.444 

10.87 

1900 

1300 

0.634 

0.656 

19.80 

650 

0.446 

0.413 

18.49 

2100 

1300 

1.040 

1.022 

22.83 

650 

1.023 

1.027 

24.01 

2200 

1300 

0.936 

1.084 

23.57 

a*x  =  [ln(l  +  a^/12)]^ 
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Here,  Bj(l)  is  the  irradiance  correlation  function.  Low-frequency 

disturbances  due  to  the  colliraation  problem  were  removed  in  Program  I 

by  a  three-section  variable  mean  and  in  Program  II  by  a  20-point  high- 

21 

pass  digital  filter.  In  addition  to  various  displays  of  the  statis¬ 
tics  of  the  log-irradiance,  the  one-dimensional  scan  power  spectra 

22 

were  calculated  by  the  usual  technique  after  prefiltering  the  data. 

The  set  c  '  one-dimensional  scan  spectra  for  each  scan  direction 
(horizontal,  vertical,  and  two  diagonal)  and  at  each  range  time  was 
then  examined  for  anisotropy.  Finally,  the  results  for  each  scan  direc¬ 
tion  were  averaged  at  a  given  range  time  to  produce  10  one-dimensional 
composite  spectra  from  an  original  120  scans. 

As  an  aid  to  further  data  analysis,  the  one-dimensional  scan 
spectra  described  above  were  fitted  to  the  form 

P(f)  =  P(fQ)  exp [Ain (f/ fc)  +  B{ln(f/f0)}2],  (17) 

where  fQ  is  an  arbitrary  reference  frequency.  Clearly,  if  B  =  0,  the 
power  spectrum  is  a  simple  power  law.  Denoting  by  a  superscript  the 
dimension  of  the  spectrum,  the  relation  between  the  isotropic  three- 
dimensional  spectra  and  its  associated  one-dimensional  spectra  is 
given  by 

A  two-dimensional  power  spectrum  is  related  to  the  one-dimensional 
quantity  by 

Through  contractual  arrangements  with  Photometries }  Inc,3  a  separate 
program  was  undertaken  for  film  scanning  and  power-spectrum  analysis 
of  the  identical  data.  The  results  of  these  efforts  are  in  direct 
agreement  with  those  reported  here . 
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> 


*(2)0O 


dY  [(Y2  +  <2)^] 


^1 
2  71 


-<X> 


dY 


1 

(Y2  +  IC2)^ 


d 

d(Y2  +  <2)^ 


(19) 


In  Appendix  A  it  is  shown  that  with  the  curve  fit,  Eq.  (17),  the  above 
equation  may  be  expressed  as 


,(2>w=^£m± 


a(Q) 


2it 


dx 


I 


k=0  CAB)  k!  dtc 


2k  X2^ 

0(|A-1+20B|;  1/2),  (20) 


where 


0  =  ln(x/Ko),  a(0)  =  A0  +  B02, 


(21) 


and  B(|A  -  1  +  2  0  B | ;  1/2)  is  the  beta  function.  Equation  (20)  is  best 
reduced  further  by  numerical  methods  in  terms  of  tabulated  i|>  functions. 
Another  approach,  and  one  used  here,  is  an  asymptotic,  or  saddle  point, 
evaluation  of  the  integral  in  Eq.  (19). 

The  average  one-dimensional  scan  spectra  are  converted  to  two- 
dimensional  log-irradiance  transverse  power  spectra  by  means  of  Eq.  (20). 
Processing  the  one-dimensional  temperature  fluctuation  spectra  by  means 
of  Eqs.  (17)  and  (18),  the  resultant  isotropic  temperature  spectrum 
divided  into  the  two-dimensional  log-irradiance  scan  spectrum  gives 
the  optical  filter  function  (OFF)  as  shown  in  Eq.  (1).  The  OFF's  for 
650  m  and  1300  m  at  particular  range  times  are  shown  in  Figs.  1  and  2. 
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RESULTS  AND  DISCUSSION 


A.  Apparent  Atmospheric  Anisotropy 

In  general,  the  one-dimensional  log-irradiance  scan  spectra  in 
horizontal  and  vertical  directions  showed  marked  differences  as  a 
function  of  frequency.  A  typical  example  is  shown  in  Fig.  6.  The 
ratio  of  horizontal  to  vertical  power  is  plotted  as  a  function  of 
frequency  for  the  1400-h  (lapse),  650-  and  1300-m  data.  There  is 
clearly  a  systematic  decrease  of  the  high-frequency  power  in  the 
horizontal  spectrum  compared  with  the  vertical  spectrum.  This  result 
might  be  interpreted  as  a  horizontal  elongation,  characteristic  of 
the  irradiance  cells,  or  as  being  due  to  an  exposure  effect  wherein 
the  shortest  wavelengths  are  transported  an  appreciable  fraction  of 
their  length  by  the  average  normal  wind  component  during  the  2-msec 
exposure.  An  estimation  of  a  wind  transfer  function  for  an  ideal 
case  is  made  in  Appendix  B.  However,  Fig.  7  shows  horizontal  and 
vertical  spectra  taken  from  a  photograph  of  a  pulsed  ruby  laser^  with 
a  30-nsec  pulse  duration.  The  anisotropy  of  the  previous  figure  is 
therefore  a  function  of  exposure  time. 

B.  Average  Scan  Spectra 

For  a  given  range  time,  the  four  individual  scan  spectra  taken 
from  the  three  photographs  were  averaged  together,  thus  reducing  the 
statistical  variation.  Samples  of  these  spectra  and  associated  beam 
photographs  are  shown  in  Figs.  8,  9,  10,  11,  12,  and  13.  The  first 
photograph-spectra  pair  gives  typical  results  for  lapse  conditions 
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Figure  6.  Ratio  of  Vertical  to  Horizontal  Scan  Power  Spectra  as  a  Function  of  Spatial 
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Average  Scan  Power  Spectra  (a)  and  Associated  Typical  Photograph  (b)  for  1900h  and  1300-m  Range 


Average  Scan  Power  Spectra  (a)  and  Associated  Typical  Photograph  (b)  for  2200h  and  650-m  Range 


(1130  h) .  The  scale  o£  scintillation  is  reduced  at  the  longer  range 
as  shown  by  the  relative  enhancement  of  the  power  at  high  frequencies. 
The  large  dark,  areas  in  Fig.  10  are  due  to  a  collimation  problem.  The 
digital  filter  applied  to  the  data  rejected  the  low  frequencies  asso¬ 
ciated  with  this  effect,  and  thus  the  low-frequency  data  do  not  appear 
in  the  result.  In  these  and  all  other  spectra,  the  high-frequency 
tail  is  shown  to  curve  up;  such  an  effect  is  due  to  optical,  film,  and 
processing  noise.  The  next  set  of  spectra,  Figs.  10  and  11,  show  beam 
photographs  at  nearly  neutral  conditions.  Clearly,  the  power  has 
dropped  in  the  scintillation  spectra  at  both  ranges,  and  sufficient 
coherence  remains  in  the  received  beam  to  produce  coherence  noise  from 
various  receiver  elements.  The  difference  of  scintillation  power  at 
the  two  ranges  appears  the  reverse  of  the  photographs  taken  during  lapse 
conditions.  This  fact  and  the  observation  that  the  peak  power  at  650  m 
is  quite  small  suggest  that  during  these  conditions  650  m  was  still 
less  than  the  transition  distance  between  geometrical  and  Fresnel  opti¬ 
cal  ranges.  Accepting  this  conjecture,  it  might  be  concluded  that  the 
inner  scale  of  turbulence  was  in  the  range  of  2.1  <  1Q  <  5.8  cm. 

During  inversion  conditions,  the  relations  of  the  two  spectra, 

Figs.  12  and  13,  are  similar  to  those  of  the  first  pair  taken  under 
lapse  conditions.  Unfortunately,  the  large  cells  in  Fig.  13  are  below 
the  cutoff  in  the  power  spectra  and  are  not  included. 

C.  Saturation  of  Statistics 

Figure  3  gives  the  measured  standard  error  of  the  log-irradiance 
versus  the  PST  value  plotted  on  the  abscissa.  The  horizontal  bars 
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represent  the  spread  in  Cn  measurements  made  by  various  techniques 
as  explained  abo%re,  while  the  vertical  bars  give  the  spread  in  data 
points.  The  data  lie  above  the  deWolf  curve  (PMST)  but  cross  it  for 
large  values  of  the  abscissa.  This  decay  of  the  measured  standard 
error  is  considered  real  and  experimentally  reproducible.  Such  be¬ 
havior  is  not  predicted  by  the  deWolf  theory,  but  may  be  due  to  finite 
beam  effects. 

D.  Optical  Filter  Function 

As  discussed  in  the  body  of  this  paper,  the  optical  filter  function 
is  derivable  from  the  data.  The  results  are  shown  in  Figs.  1  and  2. 

The  curves  marked  PST  and  BST  are  plots  of  Eqs.  (3)  and  (5)  where 
parameters  b  and  bj  were  fitted  to  our  experiment.  The  dashed  lines 
represent  measured  optical  filter  functions  in  the  frequency  range 
from  87  to  200  c/m  for  various  runs.  Data  for  which  scintillation  was 
comparable  to  system  noise  were  eliminated.  The  direction  of  the 
dashed  lines  is  considered  significant  in  distinguishing  between  the 
PST  and  the  BST  treatments.  The  frequency  range  was  chosen  so  that 
the  data  represent  the  inertial  subrange  of  turbulence  under  most 
experimental  conditions. 


CONCLUSIONS 

The  photographic  method  of  scintillation  analysis  represents 
several  distinct  advantages  over  the  prevailing  temporal  sampling 
techniques.  As  discussed  earlier,  if  a  proper  exposure  time  is  chosen 
to  freeze  the  smallest  wavelengths ,  the  data  are  sensibly  independent 
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of  the  normal  wind  speed.  Therefore,  the  uncertainty  in  the  measure¬ 
ment  of  an  average  normal  wind  speed  is  eliminated.  The  "frozen-in" 
hypothesis  need  not  he  assumed.  Finally,  the  statistical  isotropy  of 
the  beam  cross  section  can  be  investigated  readily.  Our  experiment 
shows  that  within  the  frequency  range  examined  the  refractive-index 
fluctuations  are  essentially  isotropic  even  at  a  mean  height  of  1^  m 
above  the  ground.  Clearly,  as  lower  spatial  frequencies  are  resolved, 
this  assertion  must  break  down. 

The  decline  of  the  saturation  curve  in  Fig.  3  is  considered  real. 
Such  behavior  is  not  predicted  by  the  PMST  of  deWolf  and  thus  may 
represent  a  partial  failure  of  the  initial  plane-wave  assumption  as  a 
representation  of  a  real  beam. 

Finally,  the  scan  spectra,  containing  much  more  statistical  infor¬ 
mation  than  the  standard  errors,  may  be  used  to  force  decisions  between 
competing  theories  of  light  propagation  by  means  of  the  optical  filter 
calculation.  Moreover,  the  optical  filter  function  is  relatively 
insensitive  to  the  details  of  the  index  of  refraction  spectrum  and 
consequently  is  a  universal  function  dependent  only  on  such  parameters 
as  optical  wave  number,  path  length,  and  beam  divergence.  Consequently, 
the  optical  filter  calculations  shown  as  dashed  lines  suggest  better 
agreement  with  beam-wave  theory  than  with  the  results  of  plane-wave 
predictions. 
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APPENDIX  A 


Consider  the  following  change  of  variables  in  Eq.  (19): 
X  =  K  1 (y2  +  K2) 


(A-l) 


Thus,  Eq.  (19)  becomes 


KTT 


~  exp[A(0  +  x)  +  B(0  +  x)2]  (A-2) 


(x2  -  I)'*  dx 


where 


x  =  In  x  and  e  =  ln(K/Ko). 


(A-3) 


The  derivative  with  respect  to  x  can  be  moved  out  of  the  integral  since 
^  exp[A(6  +  x)  +  B (0  +  x)2]  =  x  1  exp[A(0  +  x)  +  B(9  +  x)2]  (A-4) 
and  the  expression  is  further  simplified  with  the  chain  rule, 


d _ dK  d _ d  d 

d0  d6  d<  d<  dx 


(A-5) 


and  by  defining 


Equation  (A-2)  then  becomes 


A  =  A  -  1. 


^2\k)  =  -  iL  exp(AQ  +  B02) 

Ti 


dx 


—  exp (Ax  +  2Bx9  +  Bx2) 
(x2  -  1)* 


1 


(A-6) 

The  exponent  in  the  integrand  can  be  expressed  as  an  infinite  series 

exp  (Ax  +  2Bx0  +  Bx2)  =  exp[(A  +  2B0)x]  T 

k=0 


oo  2k  ,2k 

=  l  — ~  ~ —  2——  exp[(A  +  2B0)x]  (A- 7) 

k=0  (4B)K  k!  die 
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Thus,  Eq.  (A-fc)  becomes 


The  beta  function  exists  providing 

(A  +  2B0)  <  -1. 

Insertion  of  Eq.  (A- 10)  into  Eq.  (A-8)  gives  Eq.  (20)  of  the  text. 

It  is  also  possible  to  treat  Eq.  (A-2)  by  an  asymptotic  method. 
Let 

z  -  ( 1 B i )Ji  x, 

then  Eq.  (A-9)  becomes 


I 
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X. 


G(^)  = 


CiBl)*5 


J 

J 

0 


_ cxpfM(e)  z/dBh^-z2!  Jz 

(1  -  exp  £-2z/( 


where 


M(6)  =  A  +  2Bh. 


(A— 11) 


(A-12) 


Equation  (A-ll)  can  be  expressed  in  terms  of  a  power  series  in  z^. 
After  expansion,  Eq.  (A-ll)  reduces  to 


f 


i  -  i 

'  1  A  -h- 

x  (1  -  x)  ‘  dx 


(A-13) 


.  1  7  (-Dm  lBim 

2Jl  »•  22n 


»xp[2m  In  y  -  !■  y  -  ln(l  -  e  y)]dy. 

2  ^ 


The  first  term  in  Eq.  (A-13)  is  a  beta  function  and  the  remair.ing  series 
can  be  evaluated  by  the  saddle  point  technique.  The  integral  is  abso- 
lurely  convergent  though  the  sum  may  only  be  asymptotically  convergent. 
The  integrand  is  peaked  at  a  root  of 


£.(y)  . is.iHiaii - 1 - o. 


2  2(1  -  e'y) 


(A-14) 


An  approximate  root  of  this  equation  is  readily  extracted, 


2  (2m  -  'j) 

jH(e) ! 


providing  yQ  is  sufficiently  small  that 


1  -  e‘y  - 


(A-15) 


Conversely,  if  yQ  is  large,  then 
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(A-16) 


Am 


yo  '  |M(,  )|  +  1  * 

Both  estimates  reduce  to  the  same  value  if  2m  >  ^  and  |M(6) |  >  1. 


At  the  saddle  point, 


f (y)  :  (2m  -  -)  In  y0 


_  lM(e)[ 


2  yo 


(A-17) 


Hence 


f(y)  ,  f(y0) 

dy  e  -  e 


exp[-  f  I  f"(y0>l (y  ~  y0^dy] 


~  e 


f(yo)  ^ 


1  +  erf  {(f"(y0)y0/2)Js}]  . 


^•"(y0)l 

A  saddle  point  exists  provided  f"(y  )  <  0.  For  the  above  case, 


f"(y0>l  - 


ImB)I 

4(2ra  -  %) 


(A-18) 


Thus  an  approximation  for  G(6),  Eq.  (A-13) ,  may  be  written  as 
G(6)  =  |6(|M(0)|/2;  -|) 


(A-19) 


where  the  summation  cutoff,  w,  represents  the  convergence  "radius"  for 
the  series.  The  approximate  value  of  w  is 

.  2 


4  M(0 ) 

w=  "1 TT 


(A-20) 


The  expression  derived  for  G(0),  Eq.  (A-19),  may  be  inserted  into 

(2) 

Eq.  (A-8)  to  give  an  expression  for  4>  (k)  suitable  for  computer 

evaluation: 
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4<2>co  = 


fffal  /L\ 

2v<o  \KoJ 


[A  -  1  +  B  In (< /k  )] 


[A  +  2B  ln(K/< 


where 


Q  =  6(jM(0)J/2j 


and. 


1  w 
+  ~  [ 
/T  n^ 


cn=  I  b  |  (-1)  —Lte >  +  H)  {1  +  erf 

r(n+l)  |H(8)I  h  4  ’ 


Dn  -  |B|"(-1)"+1  - LOSt  +  3/2). 


r(n+l)  |M(0) J 


71 {1  +  erf  -  i)2]} 


qr~ “i  = 


2  )  -  ~2  j  ~  +  |M(9)[]}, 


The  <J>  functions  are  the  tabulated  derivatives  of  the  log  gamma 
In  the  event  that  B  =  0, 


Q  =  6(lM(e)  | /2;  -i) 


and  Eq.  (21)  becomes 


4>(2)(K)  =  - 


2u  r 


A-l 


A  6 (|A  -  l|/2; 


But  this  is  just  the  expression  for  a  simple  power  law  fit  of 
dimensional  scan  spectrum  with  slope  A. 


)]  Q  (A-21) 

CA— 22) 

(A-23) 

(A-24) 

(A— 25) 
functions. 

(A-26) 

(A-27) 

i  one- 
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APPENDIX  B 


It  is  useful  to  estimate  the  effect  wind  transport  of  turbulent 

has  on  the  high  spatial  frequency  components  of  the  scintillation 

spectra.  A  convenient  way  of  expressing  such  a  calculation  is  by 

means  of  an  exposure  transfer  function.  The  relation  between  the 

2 

structure  function,  D^,  of  a  quantity  f(x,y,z)  and  its  associated 
power  spectrum  in  the  constant  z  plane  is 


Df  (x-x'  ,y-y')  =  2  I  dKjd^U  -  cos[<j(x-x')  +  <2  (y-y' )]  }Ff  (*j  ,0) . 


( B—  1) 


We  will  suppose  that  a  wind,  with  velocity  vx,  blows  normally  across 
the  beam  path  and  transports  turbulence  "frozen  in"  the  wind  stream 
across  the  beam  path.  This  means,  essentially ,  that  the  current  values 
of  x,x'  (for  a  constant  structure  function)  are  related  to  the  previous 
ones  by 

x  =  x"  -  v  t;  x'  =  x"'  -  v  t'.  (B— 2) 

X  X 


Since  we  are  assuming  stationary  processes,  it  follows  that  can  only 
depend  upon  the  difference  in  the  times.  Thus, 


Df[x-x'  -  v  U-t'),y-y'] 


( B— 3) 


=  2 


dx 


x d<2 { 1  -  cosl«1(x"-x,m)  +  <2  (y-y' )  -  vx<j ( t-t ’ )] }Ff (*j ,k? ,0) 


Let  the  function  f  represent  the  instantaneous  log-amplitude,  and  let 
the  exposure  averaged  function  be  defined  as 


51 

PBEOcHlIS  PASE  BLANK 


dtdt'Df[x"-xI"  -  v(t-t') , y— y * ]  (B-4) 

where  t  is  the  exposure  time.  A  reasonable  assumption  is  that  the 
exposure  time  average  is  interchangeable  with  the  integration  in 
Eq.  (B-3) .  Consequently,  we  find 


Df (x"-x"' ,y-y' ) 


av 


j. 

~T 


(B-5) 


=  2 


Df(x"-x'",y-y')av 
T 

dK1d<2{l  -  ”7  J  dtdt '  costxj  (x"-x"')  +  K^Cy-y')  -  v^jCt-t')]} 

x  F£ (k^ , <2 » 0) • 


The  integration  over  time  is  best  accomplished  by  changing  variables 

to  the  average  and  difference  in  time.  Then  let  the  spectrum  of 

D  (x^x111  ,y-y' )  be  defined  as 
x'  J  av 


Df (x"-x"' ,y-y' ) 


av 


{1  -  cos[ic1(x"-x"’)  +  <2(y-y')]  }FfaV(icj  ,ic2>0)  dKjdic2  (B-6) 


A  comparison  of  Eq.  (B-6)  and  the  integrated  form  of  Eq.  (B-5)  gives 
the  relation  between  the  spectra 

FfaV(jc,0)  =  ^.^2t2  U  “  cos(j<’Vt)]  Ff(jc,0).  (B-7) 


The  exposure  transfer  function, 


T(k,0)  = 


77—2-2  [1  -  cos  (jc  •  vx )  ] , 


(B— 8) 


is  then  readily  identified.  The  spatial  frequency  for  which  the  trans¬ 
fer  function  has  its  first  zero  is 

f  »  500/V  cycles/m  (B— 9) 

for  the  experimental  exposure  of  2  msec.  A  fit  to  the  data  of  Fig.  6 
at  a  spatial  frequency  of  250  cycles/m  would  require  i  computed  normal 
wind  speed  of  approximately  2  m/sec,  compared  with  the  measured  wind 
speed  of  2.16  m/sec. 
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along  the  optical  path  in  order  to  estimate  the  homogeneity  of  meteorological  con¬ 
ditions.  Near  the  path  terminus,  measurements  of  wind  shear  and  temperature  lapse 
wore  taken.  In  addition,  high-speed  thermometry  techniques  were  utilized  to  compute 
onc-dimensional  temperature  spectra  as  well  as  the  thermal  structure  coefficient  C.j,. 
Data  were  gathered  during  tempera ture-lopse,  neutral,  and  inversion  conditions- 
Log- irradiance  scans  derived  from  the«y*i cal  data  were  used  to  compute  log-irradiance 
power  spectra,  variance,  and  other  staoistical  quantities.  Using  these  optical  and 
meteorological  data,  optical  filter  functions  were  then  calculated  for  spatial  fre¬ 
quencies  above  87  c/m  and  are  used  to  compare  with  current  theories.  The  saturation 
of  the  log- irradiance  data  is  again  observed,  and  the  isotropy  of  the  irradiance 
fluctuations  is  examined. 
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